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Bekenstein's Tensor- Vector-Scalar (TeVeS) theory has had considerable success as a relativistic 
theory of Modified Newtonian Dynamics (MoND). However, recent work suggests that the dynamics 
of the theory are fundamentally flawed and numerous authors have subsequently begun to consider 
a generalization of TeVeS where the vector field is given by an Einstein-TEther action. Herein, I 
develop strong-field solutions of the generalized TeVeS theory, in particular exploring neutron stars 
as well as neutral and charged black holes. I find that the solutions are identical to the neutron 
star and black hole solutions of the original TeVeS theory, given a mapping between the parameters 
of the two theories, and hence provide constraints on these values of the coupling constants. I 
discuss the consequences of these results in detail including the stability of such spacetimes as well 
as generalizations to more complicated geometries. 

PACS numbers: 04.50.Kd, 97.60.Jd, 04.80.Cc 



I. INTRODUCTION 

Modified Newtonian Dynamics (MoND) [l[ has en- 
joyed significant success as a phenomenological model of 
gravity without the requirement for dark matter (for a 
review see Q). However, the discovery of a complete 
relativistic theory generalizing MoND has proven ex- 
tremely difficult. Recently, Bekenstein's Tensor- Vector- 
Scalar (TeVeS) theory [| has received considerable at- 
tention; in the weak acceleration limit TeVeS has been 
shown to reproduce MoND, and in the Newtonian limit 
the parametrized post-Newtonian coefficients, (3 and 7, 
are consistent within limits imposed by solar system ex- 
periments (for a comprehensive review of TeVeS see [|j]). 

Despite its success, there exists mounting evidence that 
TeVeS suffers from dynamical problems. Firstly, Seifert 
Q showed that the Schwarzschild-TeVeS solution @ is 
unstable to linear perturbations for experimentally and 
phenomenologically valid values of the various coupling 
parameters. By analysing a broad variety of dynamical 
situations, Contaldi et al. Q then showed that the vec- 
tor field is prone to the formation of caustics analogously 
to Einstein- iEther theories where the vector field is given 
by a Maxwellian action 0]. Finally, Sagi 0] has shown 
that the vector field is constrained by the cosmological 
value of the scalar field in such a way that it prevents the 
scalar field from evolving and also forces it to be nega- 
tive, thereby inducing superluminal propagation of scalar 
waves. Contaldi et al. Q and Skordis [lj| have provided 
a generalization of TeVeS, whereby the vector field action 
is not given by the Maxwellian form of TeVeS, but a more 
general Einstein- jEther form. In this way the vector field 
action is the "... most general diffeomorphism invariant 
action, which is quadratic in derivatives and consistent 
with the A 2 = — 1 constraint" 0] , whereas the tensor and 
scalar field actions are identical to those of TeVeS. This 
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type of generalization to the vector field action is moti- 
vated from Einstein- iEther theories, where it was shown 
to stabilize dynamical problems. Indeed, this general- 
ization of TeVeS has been shown to resolve the caustic 
formation problem present in the original TeVeS theory 
0. 

Calculations of various quantities in generalized TeVeS 
are still in their infancy. Skordis [llj presented the cos- 
mological equations for the theory and showed that they 
are identical with the TeVeS equations up to a rescaling 
of Hubble's constant. Sagi || has looked at the PPN 
parameters within this theory and found them to not 
conflict with solar system experiments for suitable values 
of the various parameters of the theory. I continue inves- 
tigations by exploring strong-field solutions, finding the 
spherically symmetric, static solutions of the field equa- 
tions for neutron stars, neutral and charged black holes 
(sections Mil and IIV|) . Interestingly, I find all these so- 
lutions to be the same as for the original TeVeS theory 
given a simple substitution of the vector field coupling 
parameters. This implies that the work of Giannios 6] 
for neutral black holes in TeVeS, Sagi and Bekenstein 
fill ] for charged black holes and Lasky et al. [l2j for neu- 
tron stars is also correct for the generalized TeVeS the- 
ory, given the aforementioned substitution of parameters. 
Spherically symmetric neutron stars and black holes are 
therefore observationally indistinguishable in TeVeS and 
generalized TeVeS. In section fVl I discuss in detail fur- 
ther implications of this work, including generalizations 
to more complex spacetimes and stability arguments. 



II. GENERALIZED TEVES FIELD EQUATIONS 

Generalized TeVeS @, EH, as with the original TeVeS 
theory is built upon three gravitational fields; an Ein- 
stein tensor field, g^v, a normalized timelikc vector field, 
A M , and a dynamical scalar field, <p. These three objects 
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combine in the following way to give the physical metric; multiplier, A; 



9^ 



-2tp 



(g^ + A^-e^A^Av, (1) 



which is the metric measured by clocks and rods. 
Throughout this article, quantities with and without 
tildes are measured in the physical and Einstein frames 
respectively. The original and generalized versions of 
the theory differ in the vector field action. Bekenstein's 
version had a Maxwellian action, and subsequently one 
vector field coupling constant, /C. The generalized ver- 
sion has an Einstein- iEther vector field action with four 
associated coupling constants, K, K + , K 2 and K4 1 . 
Bekenstein's original TeVeS theory is recovered by let- 
ting K + = K 2 = K4, = and K = K. 

The modified Einstein field equations can be expressed 
analogously to TeVeS as 
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Here, is the Einstein tensor associated with the Ein- 
tein frame and T^„ is the physical stress-energy tensor. 
The tensor, t^ u , has the same definition as that in the 
original theory 
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where k is the scalar field coupling constant, £ is a positive 
scalar of dimension length, /j, is a function associated with 
the MoND acceleration scale, and T (/x) is a function 
that is not predicted by the theory. Only the strong-field 
limit of the theory is considered in this article, whereby 
an excellent approximation is [/, = 1 (for more details 
see [1, H, 0, EH)- Moreover, in this case one finds that 
T diverges logarithimcally, but is exactly cancelled out 
in the field equations implying one is free to ignore this 
function Q- It is worth noting that T only diverges 
in the strong-field limit considered herein, implying this 
divergence is not a general concern for the theory. 

As mentioned, the extra degrees of freedom attributed 
to the generalized version of TeVeS appear in the vector- 
field action, which manifests itself partly in the M „ term 
of equation ((2|) (for details of these derivations see [7, 
03). For clarity in the expressions, can be expressed 
as a linear sum of terms associated with each coupling 
parameter, plus one term associated with the Lagrange 



1 I follow the notation of Sagi which is different to that of 
Skordis [lOl who denotes the four parameters by Kg, K+, Kq 
and Kj\. These are related as follows; K = Kg = (ci — c^)/2, 
K + = K + = (cj + c 3 )/2, K 2 = K = c 2 and K 4 = K A = -c 4 , 
where the Ci are those of the Einstein- jEther theory (see Jacobson 

El)- 



' JJ.U 



where 



--K I F afl F a u - -g^F a(3 F 



a/3 



Q Q a _ 1 rr ca/3 



Q% := K 2 



S^V Q {A a Vf}AP) - A {i y u) (V a A a ) 



^g^ a A a V p AP 



Ap,A u + A a A^Vy) A a 



V Q ^A A^Ay 



1 



9iiuA a A a 



e 



XA^Ay. 



(4) 
(5) 

(6) 

(7) 

(8) 
(9) 



Here, F^ := V [v A^, S := V( U A^ and A^ := A Q V Q A M . 
Additional to the modified Einstein equations, one has 
the vector field equation given by 

KW a F^ a + K + V a S a » + K 2 V (V a A a ) + XA" 
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which again reduces to the original TeVeS vector field 
equation when K + = K 2 = K 4 = and K — K,. Finally, 
the scalar field equation is 

V [/x {g af} - A a AP) V a ip] 

= kG [g ap + (1 + e" 4 *) A a A lS ] f af} , (11) 

which is equivalent to that of the original TeVeS theory. 



III. NEUTRON STARS AND NEUTRAL BLACK 
HOLES 

Giannios [6(] first considered static, spherically sym- 
metric vacuum solutions of the original TeVeS field equa- 
tions. He found two branches of solutions; one where the 
vector field has only a non-zero timelike component, and 
one where the vector field has both non-zero radial and 
temporal components. Lasky et al. [HJ then considered 
spherically symmetric neutron star solutions of the origi- 
nal TeVeS field equations where the vector field only has 
a non-zero temporal component. The black hole solution 
of Giannios 3 can be recovered from the neutron star 
solutions of Lasky et al. H3| by simply letting the matter 
variables vanish throughout the spacetime. I therefore 
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begin here by deriving the neutron star solutions and sub- 
sequently the black hole solutions. Moreover, throughout 
this section I consider a vector field whichpoints in the 
temporal direction, for which Bekenstein [3[ has shown 
support for the general case of static spacetimes (for a 
discussion of generalizing this work to a more general 
vector field ansatz see section |V|) . 

A spherically symmetric, static spacetime can be ex- 
pressed in the Einstein frame as 
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where dfl 2 := d9 2 + sin 2 8d<fi. As mentioned, the vector 
field is chosen such that is has only a non-zero temporal 
component. Normalization of this implies 

= 5^ t e- u/2 , (13) 

further implying that the physical line element is 

g a pdx a dx 13 = -e" +2v dt 2 + e' 2v (e^dr 2 + rW) , (14) 

where symmetries also imply tp — tp(r). 

Stress-energy tensors in TeVeS and its generalization 
are defined in terms of the physical frame. Therefore, a 
perfect fluid takes the form 
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where 



p and P are respectively the four-velocity, 



energy-density and pressure of the fluid as measured in 
the physical frame. Conservation of the stress-energy 
tensor is given by VqT"^ = 0, where V M is the unique 
metric connection associated with the physical metric. 
The radial component of this equation implies the pres- 
sure gradient is expressed as 



(16) 



where a prime denotes differentiation with respect to the 
radial coordinate. 

The scalar field equation (fTTj) can be expressed in terms 
of the metric coefficients, and once integrated to give 
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where the scalar mass function, M v {r), has been defined 
as 

M v (r) := 4tt J (p + 3PJ c v ' 2+ ^ 2 - 2 ^f 2 df. (18) 

Given the symmetries of the problem and the vector field 
ansatz, the radial component of the vector field equa- 
tion is now trivially satisfied, and the temporal compo- 
nent provides an equation for the Lagrange multiplier, A. 
Therefore, this is essentially a periphery equation as the 



Lagrange multiplier is used in the modified field equa- 
tions in the 0*, term of equation The only remain- 
ing equations are the modified Einstein field equations 
@. After much work, the tt, rr and 86 components can 
be respectively shown to be 



rC + e c - 1 = 8TrGpe c - 2v r 2 + — (K + K+ - K A ) 



v' 2 „ ,/ 4i/' 

\ I — v'C + 2" + 

2 r 



47rH 



-tp- 



(19) 



rv 1 + c c + 1 =8irGe c - 2 ' fi Pr 2 — (K + K + — K 4 ) —v' 2 

8 

47rr 2 



ii 



(20) 



- (2v' - 2(' - ru'C + 2riv" + rv 12 ) = %TiG^- 2ip Pr 2 

+ !-(K + K + -K 4 )is /2 -^-p' 2 . (21) 



k 



A striking similarity between this system of equations 
and the equivalent system in the original TeVeS theory 
presented in Lasky et al. [Hj], their equations (21-23), is 
now observed. In fact, the system of equations presented 
for the original TeVeS theory is exactly the same, where 



K + K+ - K 4 = K, 
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in the above equations. That is, given the neutron star 
solution of the original TeVeS theory, one derives the 
solution of generalized TeVeS by making the substitu- 
tion (|2"2"|) . The spacetime of neutron stars in both theo- 
ries is therefore identical. Moreover, allowing the matter 
terms to vanish in equations (|16ti21[) . and therefore re- 
covering the equations governing the generalized TeVeS 
Schwarzschild spacetime, one can show that the structure 
of the spacetime is exactly equivalent to that presented in 
Giannios @ for the original TeVeS theory. To show this, 
one must first perform a coordinate transformation given 
by r — r(R) such that r = e^ R ^ 2 R to put the system 
in the same set of coordinates used by Giannios. Then, 
allowing the matter terms to vanish, the system of equa- 
tions becomes exactly equivalent to equations (28-30) of 
Giannios [6j, given the expression in (f2"2"|) . That is, the 
structure of black hole spacetimes, under the symmetry 
assumptions posed hitherto, is exactly the same as that 
of the orignial TeVeS theory. 

There is no need to further derive solutions of these sys- 
tems of equations for both neutron stars and black holes, 
as the above result implies the results of Lasky et al. [l2j 
and Giannios [6j hold for the generalized TeVeS theory, 
providing equation (|22[) is utilized. One observes from 
equations (|16H21[) that the individual coupling constants 
do not appear independently of one another, but rather 
neutron stars and black holes can only be used to con- 
strain the combination of parameters, K + K + — K 4 . As 
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with Lasky et al. [12J, we can use observations of neutron 
star masses to constrain these coupling constants. In that 
paper, we gave a conservative estimate of K. < 1 based 
on the existence of neutron stars with masses of at least 
~ 1.5M Q . Therefore, using the same conservative esti- 
mate the combination of the vector coupling constants is 
constrained to K + K + — K4 < 1. 

In some ways the relationship between the two TeVeS 
theories given in (|22l) is to be expected. As was discussed 
in section |TT1 one recovers the original TeVeS theory by 
setting K + = K2 = K4 = and K = K, which is an 
obvious subset of solutions offered by the above equa- 
tion. However, (|22|) does not imply that this has to be 
the case, but rather any combination of K, K + and 
may hold such that (|22[) is true. Moreover, these space- 
times are completely independent of the parameter Ki. 
I discuss the consequences of these result in considerably 
more detail in section fVl 
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Moreover, the scalar field equation (fTT|) can be shown to 
reduce to 



IV. CHARGED BLACK HOLES 

Sagi and Bekenstein [1 ll ] considered spherically sym- 
metric, charged black holes in the original TeVeS theory, 
where the vector field is again aligned with the temporal 
direction. In this section I reproduce those calculations 
in the context of the generalized TeVeS theory. 

The ansatz and the symmetries applied herein imply 
the metric and vector field are given by equations (|12j) 
and (fT3]) respectively. The stress-energy tensor is now 
taken to be that of an Einstein-Maxwell field as measured 
in the physical frame; 
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where is the electromagnetic Faraday tensor in the 
physical frame. Conservation of stress-energy then leads 
to the Maxwell equations, V tt P° = 0. The symmetries 
of the problem imply the only non-zero contribution to 
the Faraday tensor is the T r t — —Tt r component [TT[ | . 
Moreover, only the temporal component of Maxwell's 
equations is non-zero, which is once integrated to give 
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Here, Q is a constant of integration which is the charge 
of the black hole. As with section Hill the only non-zero 
contribution from the vector field equation (|10p is from 
the temporal component, which gives the Lagrange mul- 
tiplier expressed as a function of the metric coefficient. 
Evaluating the modified Einstein equations ^ and sub- 
stituting the Lagrange multiplier back through gives the 
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Sagi and Bekenstein [Tl| derived the same set of equa- 
tions in the original TeVeS theory, although they used 
isotropic coordinates as opposed to Schwarzschild coor- 
dinates. Performing the same coordinate transformation 
as section IIII1 i.e. r = r(R) such that r = e^^ R ^ 2 R 
where the Einstein metric has component g^ji = e^ R ' , 
and also using equation (j22[) . one finds that equations 
(|25H28[) given above are identical to equations (49-52) of 
Sagi and Bekenstein That is, given the charged 

black hole solutions of Sagi and Bekenstein [ll| , one can 
derive the solution of generalized TeVeS simply making 
the substitution (|22[) , implying the spacetimes in the two 
theories are identical. 



V. IMPLICATIONS AND GENERALIZATIONS 

Herein I have shown that solutions of the field equa- 
tions within Bekenstein's original TeVeS framework for 
spherically symmetric neutron star s 11211 . neutral black 
holes @ and charged black holes [ill] are exactly the 
same solutions for the generalized TeVeS theory [3, [l(| , 
given the substitution (]22[) . An immediate question one 
must ask is; how much can these results be generalized? 
That is, if more complex solutions of the original TeVeS 
field equations are found, will the same solution of the 
generalized TeVeS equations exist given the substitution 
(|22"f ? The answer to this question is a simple "no" . In 
fact, a quick calculation reveals that this is not even the 
case in more complex spherically symmetric spacetimes. 
For example, consider the same metric ansatz utilized 
in this paper, and allow the vector field to have a non- 
zero radial and temporal component. One can then show 
that the divergence of the vector field is non-zero, i.e. 
V a j4 Q ^ 0, and consequently that the Ki vector field 
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coupling now has a non-zero contribution to the modified 
Einstein equations. Therefore, the simple mapping given 
in (|22p can no longer hold between the original TeVeS 
theory and the generalized version. Moreover, one can 
immediately see that when time dependence or rotation 
is included into the spacetime, the K2 term will gener- 
ally have non-trivial contributions due to the non-zero 
divergence of the vector field, and subsequently solutions 
of the original TeVeS theory will differ markedly from 
those of the generalized TeVeS theory. 

Implications for observations of neutron stars and 
black holes given the above information is abundantly 
clear. Indeed, this article implies that spherically sym- 
metric neutron stars and black holes in the general- 
ized TeVeS theory, where the vector field is temporally 
aligned, are indistinguishable from those of the original 
TeVeS theory. Meanwhile, the above discussion suggests 
that observations of rotating black holes and neutron 
stars will differ between the theories, although when slow 
rotation approximations are valid, one would expect the 
observations to not differ significantly. 

The relation provided in (|22|) is not unique to this ar- 
ticle. In fact, Skordis showed that the combination 
K + K + — K4 plays the same role in cosmological per- 
turbations of the generalized TeVeS theory as /C in the 
original theory. Moreover, he showed that to keep the en- 
ergy density of cosmological perturbations positive, the 
constant must satisfy < K + K + — K4 < 2, which 
is consistent with phenomenological constraints placed 
of < K, < 2 in the original TeVeS theory p], El. 
This is further consistent with Einstein-jEther theories, 
where this combination of constants (commonly denoted 
C14 = ci + C4) must also satisfy C14 < 2 such that New- 
ton's constant is positive. 

An interesting question regarding the stability of the 
generalized TeVeS theory now arises. Indeed, the first 
blow was dealt to the stability of the original TeVeS the- 
ory when Seifert [5J found that the Schwarzschild-TeVeS 
solution 6] is unstable. As the structure of the spacetime 
in the generalized theory is the same, does this necessarily 
imply that the black holes are also unstable? It is possi- 
ble that this is not the case for the following reason; in- 
stabilities are a consequence of dynamical perturbations 
of the spacetime. Therefore, while the overall structure 
of the background spacetime is the same for the original 
and generalized TeVeS theories, the field equations differ 
due to the extra contributions of the vector field action, 
and therefore the dynamics of the perturbations will also 
differ. Whilst full calculations of this are still required, 
one can draw on experience from Einstein- jEther theory. 
There it has been shown that the vector field develops 
caustic singularities as it falls into the potential well of a 



gravitational field when the vector action is Maxwellian 
Q]. For this reason, generalizations of Maxwelllian vector 
fields are generally considered in modern Einstein- iEthcr 
literature. Indeed this was the original motivation for 
generalizing TeVeS in such a way, although I emphasize 
the need for these arguments to be verified in the gener- 
alized TeVeS theory by way of rigorous calculations. 

Whilst it is true that general perturbations will behave 
differently due to the different field equations of the the- 
ory, a certain subset of perturbations will behave in the 
same manner. For example, Sotani [14j | studied perturba- 
tions of neutron star solutions in the original TeVeS the- 
ory by assuming a Cowling approximation. Here, the per- 
turbation equations are simplified considerably by only 
dealing with perturbations of the fluid variables, and not 
the background spacetime (i.e. the metric, vector field 
and scalar field remain fixed). As the fluid equations 
in the physical frame are identical for the original and 
generalized TeVeS theories, the results of Sotani [l4| also 
hold true in generalized TeVeS, where again the mapping 
between the parameters is given by (|22"|) . 

In summary, in this article the equations govern- 
ing spherically symmetric, static spacetimes represent- 
ing neutron stars, neutral and charged black holes in the 
generalized TeVeS theory have been presented. These 
turn out to be equivalent to those of Bekenstein's origi- 
nal TeVeS theory H given by Lasky et al. [12] , Giannios 
[6| and Sagi and Bekenstein , where the coupling con- 
stants of the two theories are related by (|22| . Therefore, 
the solutions of those equations and hence the structure 
of the spacetimes presented by those authors is the same 
as for the generalized theory, and indeed these objects are 
observationally indistinguishable in this regime. I have 
used this to show that the combination of vector coupling 
constants can be constrained to K + K + — K4 < 1 based 
on conservative neutron star observations. Moreover, 
whilst problems with the original TeVeS theory were ini- 
tiated by studies of spherically symmetric black holes, it 
is not expected that this will be the case for the gener- 
alized theory, given that perturbations and dynamics of 
such spacetimes will behave differently. 
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